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ABSTRACT 



We report about the search for short-term variability in the high-energy y-ray energy band of three flat-spectrum radio quasars 
(3C 454.3, 3C 273, PKS B 1222+216), whose flux at £ > 100 MeV exceeded the value of 10"' ph cm^^ s"' for at least one day. 
Although, the statistics was not yet sufficient to effectively measure the characteristic time scale, it allowed us to set tight upper limits 
on the observed doubling time scale (< 2 - 3 hours) - the smallest measured to date at MeV energies ~, which can constrain the size of 
the y-ray emitting region. The results obtained in the present work favor the hypothesis that y rays are generated inside the broad-line 
region. 
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The latest generation of ground-based Cerenkov telescope, like 
HESS and MAGIC, has detected flux changes on a few minutes 
time scale in the BL Lac Objects PKS 2155-304 (Aharonian et 
al. 2007) and Mkn 501 (Albert et al. 2007). These episodes, al- 
though exceptional, have severely challenged the common as- 
sumptions about the size and location of the emitting region. 
According to these models, the emitting region of size r > r^, - 
GM/c^ (being rg the gravitational radius of the central black hole 
with mass M) is located at a distance R ~ r/tp from the central 
spacetime singularity, where i/r is the aperture of the jet (gener- 
ally i/f ~ 0.1 - 0.25, see Ghisellini & Tavecchio 2009, Dermer et 
al. 2009). Changes in the emitted radiation from this blob occur 
on a characteristic observed time t/(1 H-z) > r/c6, where S is the 
Doppler factor of the jet. 

It is worth noting that the dissipation zone should be suffi- 
ciently far from the central black hole to avoid pair production 
due to the interaction with the photons emitted mainly by the 
accretion disc and the broad-line region (BLR). BL Lac Objects 
are known to have a photon-starved environment: the accretion 
disc in these sources is weak, likely an advection-dominated 
disc (ADAF), which is also inefficient in ionizing the BLR. 
Therefore, the variability of a few minutes observed in these BL 
Lac Objects can be explained by assuming very small blobs at 
small R: the paucity of soft photons makes it possible for GeV 
photons to escape without pair production. 

This is no more possible in the case of FSRQs, where discs 
with accretion rates of a few tens of percent of the Eddington 
rate are present, which in turn efficiently photoionize the BLR, 
generating broad and strong optical emission lines. This rich en- 
vironment (disc, BLR) makes the cooling of relativistic electrons 



more efficient through the external Compton (EC) processes, 
but, on the other hand, if the emitting blob is too close to the 
accretion disc to be sufficiently small for the short variability, 
then the pair production could severely absorb the GeV photons. 

Some researchers have proposed different solutions to this 
problem, which are basically relying on the increase of the 
Doppler boosting or some special jet structure (Begelman et 
al. 2008, Finke et al. 2008, Ghiseflini & Tavecchio 2008, 
Giannios et al. 2009). The latter, based on the hypothesis that 
there are very compact relativistic blobs within the jet ("nee- 
dle/jet" model, Ghisellini & Tavecchio 2008; "jet-in-jet" model, 
Giannios et al. 2009), can be adapted to both BL Lac Objects and 
FSRQs, although it seems an ad hoc solution at the present stage 
and, perhaps, it is necessary to better assess this type of solution 
by comparing with more and more observations. Therefore, the 
most critical aspect of these theories would be the discovery of 
variability of a few minutes in the case of FSRQs. 

The shortest time scales measured in FSRQs is of about half 
an hour, as observed in the hard X-ray emission (20-40 keV) of 
NRAO 530 in 2004 February (Foschini et al. 2006). However, 
given the energy range, it would be possible that this short flare 
could have been generated by the X-ray corona, although the 
emitted power (~ 8 x 10'*'' erg s"') seems a bit large to be due 
to an unbeamed component. In addition, an increase of the radio 
polarization around the time of the flare suggests that it was in- 
deed related to the jet, although the coarse observation sequence 
did not allow to set tight constraints and doubts shadow the con- 
clusions. 

Anyway, if short time scales could be observed in FSRQs at 
hard X-rays, which are sampling low energy part of the inverse- 
Compton hump of the spectral energy distribution (SED), then 
there is the possibility that similar variability is present also in 
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the high-energy part (E > 100 MeV). At these energies, there is 
no more doubt that the emission is coming from the jet, since no 
other structures are known to generate such energetic photons at 
cosmological distances. 

The known flux variations at high-energy y rays {E > 
100 MeV) in FSRQs are generally in the range of a few hours, 
in agreement with the current paradigm. The best examples 
measured by CGRO/EGRET during nineties are PKS 1622-29 
(~ 4 hours, Mattox et al. 1997) and 3C 279 (~ 8 hours, 
Wherle et al. 1998). The early results from Fermi satellite ob- 
tained in the latest years have confirmed such time scales: 
PKS 1454-354 (Abdo et al. 2009), PKS 1502-H106 (Abdo et 
al. 2010a), PKS B 15 10-089 (Tavecchio et al. 2010), 3C 454.3 
(Foschini et al. 2010, Tavecchio et al. 2010, Ackermann et al. 
2010) and 3C 273 (Abdo et al. 2010b). 

The availability of the Fermi satellite provides a plethora of 
y-ray data where to search for short variability events. The Large 
Area Telescope (LAT, Atwood et al. 2009) onboard Fermi rep- 
resents the state of the art of y-ray space instruments, with an 
increase of sensitivity of a factor ~ 20- 30 with respect to its 
predecessors. Fermi operates in scanning mode, i.e. it performs 
an all-sky survey every 3 hours (two orbits). This is the first con- 
tinuous monitoring of the high-energy y-ray sky and offers the 
unique possibility to study the blazar population, the duty cycle 
of individual sources, and to catch the most powerful outbursts. 
The bad side of the thing is that the scanning mode hampers the 
probing of subhour variability, because the source is not always 
at the LAT's boresight, where the instrument has its best perfor- 
mance. A tentative to get over this obstacle out has been done in 
2010 April, when 3C 454.3 underwent an intense ourburst with 
flux above 100 MeV in excess of 10"^ ph cm"^ s"^ During this 
episode, Ferm//LAT performed a pointed observation staring at 
the FSRQ (2010 April 5-8). The latter, however, did not col- 
laborate and remained almost constant with poorly significant 
variations (~ 3 hours) and began to be variable only toward the 
end of the special observation (Foschini et al. 2010, Ackermann 
et al. 2010). 

Therefore, we decided to expand our search to other - 
possibly more active - periods and sources. We searched for 
other cases of FSRQ with high-energy y-ray flux exceeding 
10"^ ph cm"^ s"' for at least one day, in order to have the best 
available statistics. We found three FSRQs fulfilling this crite- 
rion, which were 3C 454.3 (z = 0.859), 3C 273 (z = 0.158), 
and PKS B1222H-216 (z = 0.432). We report here the results 
of this analysis. There were two more cases that could have 
been of interest, but they were discarded. PKS B 15 10-089 and 
PKS 1830-21 1 have exceeded the threshold flux for a few hours 
(Ciprini et al. 2010, Tavecchio et al. 2010), but not for at least 
one day, and therefore we did not consider them. Moreover, at 
the time of writing this work (Christmas 2010), 3C 454.3 ex- 
ceeded again the threshold, but we limited our data set at the end 
of 20 10 November. 



2. Data Analysis 

As stated above, we have searched for very high fluxes sources 
(Fy > 10"^ ph cm"^ s"' averaged over one day; E > 100 MeV). 
We have identified three candidates: 3C 454.3, which has ex- 
ceeded the flux threshold three times to date (2009 December, 
2010 April and November), PKS B1222+216 (2010 May-June) 
and 3C 273 (2009 September). Fermi/LAT data for the above 
mentioned sources and time periods were downloaded from the 
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Fig. 1. Light curves of 3C 273 (E > 100 MeV) between 2009 
September 15 and 23. Time starts on 2009 August 31 (MJD 
55074), so to have in the abscissa the days of September. Time 
bins are of the order of a few thousands of seconds and are too 
small to be visible. 

Fermi Science Support Center at HEASARCQ. The selected data 
were screened, filtered, and analyzed as described in greater de- 
tail in Foschini et al. (2010), but with a more recent version of 
the LAT Science Tools (v. 9. 18.6) and the corresponding 
calibration files. 

Moreover, while in Foschini et al. (2010) we could build time 
bins smaller than the good-time intervals (GTI) because of the 
better LAT performance during pointed observations (the col- 
lected counts are greater than in survey mode by a factor ~ 3.5, 
see also Ackermann et al. 2010), this is no more possible when 
analyzing LAT data from scanning mode observations, which in 
turn constitute the majority of the data analyzed in the present 
work. Searching for the best trade off between a small time 
bin and the need of significant statistics in each bin, we have 
found that the best option is to have time bins equal to the GTI, 
which generally are of the order of a few kiloseconds (a bit less 
than one orbit). Shorter bins could suffer of artifacts as indicated 
in the caveats listed in the LAT Science Tools web pages at 
HEASARC; larger bins would wash out the variability. 

Last, but not least, we have required that the flux in each 
bin was at least a factor 2 greater than its error, i.e. that there 
were sufficient events to correctly and significantly reconstruct 
the source flux. 

The extracted light curves are displayed in Figs. ([T]-[3]). We 
have noted that 3C 454.3 reached its peak flux of (1.0 ± 0.1) x 
10""* ph cm"^ s"' {E > 100 MeV, coiTesponding to a luminosity 
of about ~ 3 X 1 0^" erg s" ' ) on 20 1 No vember 20 between 1 :45 
and 03:03 UTC (source ontime ~ 4.7 ks). During this time, LAT 
detected 1 10 events from the blazar (22cr detection). This is the 
greatest y-ray flux ever detected to date from an AGN. It is worth 
noting that, contrary to previous observations at high-energy y 
rays, this time 3C 454.3 displayed significant spectral changes 
in the y-ray energy band (Abdo et al. 201 1). 

Once prepared, we have scanned all the light curves search- 
ing for the minimum time of doubling/halving flux: 



F(t) = F(fo) ■ l-'-'-'"^''' 



(1) 



where F(t) and F(to) are the fluxes at the time f and fo, respec- 
tively, and T is the characteristic time scale. We have required 
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Fig. 2. Light curves of PKS B1222+216 {E > 100 MeV). (left panel) Light curve in the period 2010 April 29 (MJD 55315) and 
June 20. (right panel) zoom of first ten days of the curve. Time bins are of the order of a few thousands of seconds and are too small 
to be visible. 
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Fig. 3. Light curves of 3C 454.3 (£ > 100 MeV). (top left panel) 2009 December: time starts on 2009 November 30 (MJD 55 165), 
so that the days indicated in abscissa corresponds also to the days of December (top right panel) 2010 April: time starts on 2010 
March 3 1 (MJD 55286), so that the days indicated in abscissa corresponds also to the days of April. The time region between the 
two dashed lines is that when the pointed observation was done (2010 April 5-8). It is clearly evident that the error bars are smaller 
than those in scanning mode, as a consequence of the better performances of the instrument, (bottom left panel) 2010 November: 
time starts on 2010 October 31 (MJD 55500), so that the days indicated in abscissa corresponds also to the days of November. 
(bottom right panel) 2010 November: zoom of the figure of the bottom left panel centered in the period of the highest activity (2010 



November 16-22), when 3C 454.3 reached the peak of - 
and are too small to be visible. 



10 ph cm s . Time bins are of the order of a few thousands of seconds 
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Table 1. Summary of characteristic times found. Times to and t are in [MJD]; fluxes are in units of [10^^ ph cm"^ s"']; the 
significance of the flux difference is in [cr]; the absolute values of the observed characteristic time scale |r| and the intrinsic one 
kinti - + z) are in hours. The last column indicates a rise (R) or a decay (D) time. 



Source 


Period 


to 


t 


F(to) 


Fit) 


Signif. 


\T\ 


|T,nll 


R/D 


3C 273 


2009 September 


55094.88299 


55095.02236 


0.28 + 0.14 


0.87 + 0.06 


4.2 


< 4.9 


< 4.2 


R 


1222+216 


2010 May- June 


55315.76870 


55315.83500 


0.82 + 0.18 


0.18 + 0.07 


3.5 


< 2.3 


< 1.6 


D 






55315.83500 


55315.90130 


0.18 + 0.07 


0.80 + 0.20 


3.0 


< 2.5 


< 1.7 


R 






55322.52962 


55323.59231 


0.08 + 0.04 


0.51+0.08 


5.4 


9.5 + 5.3 


6.6 + 3.7 


R 


3C 454.3 


2009 December 


55167.16269 


55167.21717 


0.85 + 0.37 


2.53 + 0.56 


3.0 


< 3.4 


< 1.8 


R 






55171.33812 


55171.40897 


2.06 + 0.37 


0.73 + 0.23 


3.6 


< 3.6 


< 1.9 


D 




2010 November 


55511.48731 


55511.55361 


1.32 + 0.28 


0.45 + 0.22 


3.1 


< 5.0 


< 2.7 


D 






55514.12842 


55514.19822 


0.57 + 0.20 


1.63 + 0.35 


3.0 


< 4.1 


< 2.2 


R 






55514.47072 


55514.53701 


1.94 + 0.39 


0.76 + 0.21 


3.0 


< 4.7 


< 2.5 


D 






55522.66673 


55522.73602 


3.42 + 0.46 


1.27 + 0.59 


3.6 


<4.3 


< 2.3 


D 






55524.28281 


55524.34911 


1.68 + 0.27 


3.71+0.62 


3.3 


< 4.8 


< 2.6 


R 



that the difference in flux is significant at 3cr level (at least) and 
we have selected only the results with |t| < 5 hours (including 
the uncertainties). The results are displayed in Table fllrl. 

In the case of 3C 454.3, we note that the less variable pe- 
riod is that of 2010 April, just when the pointed observation 
was done. This was already noted in Foschini et al. (2010) and 
Ackermann et al. (2010), although the former reported a ~ 3cr 
change in flux at the end of the pointed observation (in the early 
hours of 2010 April 8, when there was a short flare; see Fig. |3] 
top right panel) with t ~ 3 hours, which is not found in the 
present analysis, because - as expected - a greater time bin 
(1.8 ks in Foschini et al. 2010; 4.8-5.7 ks in the present work) 
smoothed the variability (the measured change in flux is now at 
~ 2.5cr level). Therefore, since the pointed observation has been 
already studied in the above cited works and the analysis per- 
formed in the present work cannot give better results, because 
it is best suited for observations in scanning mode, we have de- 
cided to drop this period and we report the data in the figures for 
the sake of completeness. 



3. Discussion 

The masses of the three sources are of the order of some X lO^M© 
(see Bonnoli et al. 201 1, Fan & Cao 2004, GhiselUni et al. 2010), 
which means gravitational radii rg ~ (0.5 - 1) x 10^"* cm and a 
minimum characteristic time r^/c ~ 0.5 - 0.9 hours. The min- 
imum size of the emitting blob of plasma as inferred from the 
upper limit of variability is r < c6Ti„t ~ (20 - 30)r,,, by con- 
sidering a typical 6 = 10 (Ghisellini et al. 2010). Moreover, 
the size r of the emitting blob is linked to its distance from 
the central engine by the relationship r ~ ifrR, where iff is 
the semi-aperture angle of the jet, generally about 0.10 - 0.25 
(Ghisellini & Tavecchio 2009, Dermer et al. 2009). The val- 
ues of r calculated from the measurements of variability in the 
present work are r ~ (1 - 3) x 10'^ cm, which then imply 
R ~ (1 -8)x 10'^ cm, smaller but comparable with the BLR size, 
where most of the dissipation occurs in FSRQs (cf Ghisellini et 
al. 2010). Therefore, although the characteristics time scales re- 
ported in the present work are the shortest ever measured in the 



^ It is worth noting that there are some flares in the curves resolved 
with a few points. In this case, it would be possible to fit these points 
and calculate the t with smaller errors. For example, in the case of 
PKS B1222+216 it is possible to measure t = 1.5 + 0.6 hours (fit of 
4 points). However, given all the uncertainties and the low statistics, we 
prefer to adopt the most conservative approach outlined above. 



MeV energy range to date, they are still satisfying the conditions 
on the spatial dimensions suggested from the standard paradigm. 

We can try to understand if these upper limits on t could 
be useful in distinguishing between the location where most of 
the dissipation occurs. This is particularly important, because 
during the latest years there was and still is a lively intriguing 
debate on this topic. On one side, there is a group of researchers 
supporting the hypothesis that the y rays in the sources powered 
by the EC process are generated by the interaction with seed 
photons from the molecular torus (infrared photons), which is 
placed at some parsecs from the central engine (e.g. Abdo et 
al. 2010c; Marscher et al. 2008, 2010; Sikora et al. 2009). On 
the other side, another group of researchers is instead favoring - 
with different arguments - the hypothesis of a dissipation inside 
the BLR, i.e. less than one tenth of parsec (e.g. Finke & Dermer 
2010, GhiselHni et al. 2010, Poutanen & Stern 2010, Tavecchio 
et al. 2010). 

In order to distinguish between the two hypotheses, we have 
followed the argument on the electrons cooling time as explained 
in Tavecchio et al. (2010). Basically, the observed cooling time 
jobs ^.^jj calculated as: 

"4crT£/'\ v°b^5 

where me is the electron rest mass, (Tt is the Thomson cross 
section, F is the bulk Lorentz factor, v"^^^^ and v°^^ are the peak 
frequencies of the seed photons (2 x lO'^ Hz for the broad- 
line region; 3 x lO'^* Hz for the molecular torus) and of the 
inverse-Compton (100 MeV= 2.4 x 10^^ Hz). U' is the en- 
ergy density of seed photons (in the comoving frame), which 
is ~ 3.8 X I0~^ X erg cm"^ for the broad-line region and 
~ 3.0 X 10 X erg cm"^ for the infrared torus (see Ghisellini 
& Tavecchio 2009). The values calculated by using Eq. ^ have 
to be compared with the characteristic times observed during the 
declining phases of the outbursts (see the D values of the last 
column in the Table [TJ. 

By assuming the typical value of F ~ 5 ~ 10 (e.g. Ghisellini 
et al. 2010), Eq. (O becomes: 

t^^^ ~ 0.67 X V(l +z) hours (3) 

jtoms ^ J ^ ^(1 +z) hours (4) 

By substituting the redshifts of the sources, the observed 
cooling times span from a bit less than one hour in the case 
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of the broad-line region and more than 10 hours for the torus. 
Therefore, the observed upper limits of a few hours (see |t| in 
Table [1]) of the characteristic time scales favor the BLR (sub- 
parsec) location of the dissipation zone. 

It would be possible to reconcile the cooling time of EC 
based on the infrared seed photons from the molecular torus if 
assuming a large bulk Lorentz factor. For example, if F ~ 5 ~ 30, 
then Eq. (|4]i becomes: 

~ 1.0 X V(l +2)hours (5) 

However, this is an unlikely hypothesis, in conflict with the 
values of the jet apparent speed as measured with high-resolution 
radio observations. Lister et al. (2009), by analyzing a large 
set of data obtained with high-resolution VLBA observations at 
15 GHz in the period 1994-2007, have reported that the greatest 
values of the apparent speed are generally measured at positions 
very close to the core, at the limit of the angular resolution of the 
instrument (~ 1 mas). In the cases of the three sources studied in 
the present work, they reported the values of jS™'' - 21, 13, 14 
for PKS B1222+216, 3C 273 and 3C 454.3, respectively, which 
corresponds to Doppler factors roughly in the range 5 ~ 9 - 15 
(calculated by assuming jSapp ~ ^26Y). The locations of the fea- 
tures used to perform these measurements are at ~ 5 - 15 pc, 
which are the typical distances of the infrared torus. A Doppler 
factor greater by more than a factor 2 seems to be unlikely, al- 
though it cannot be excluded a priori: we are studying outstand- 
ing high fluxes at y rays and, therefore, we cannot avoid thinking 
that also radio observations (if any) following these strong out- 
burst could result in exceptional values of ySapp. 

4. Final remarks 

We have searched for the shorter time scale in the high-energy 
y-ray emission from FSRQs. To have the best trade off be- 
tween the smallest time bins and sufficient statistics for a sig- 
nificant detection, we have analyzed the Ferwi/LAT data of 
three FSRQs, whose flux above 100 MeV has exceeded the 
threshold of 10"^' ph cm"^ s '. We have found three sources: 
PKS B1222+216, 3C 273, and 3C 454.3. We were able to set 
the tightest upper limit on the observed characteristic time scale 
to date, which were of the order of < 2 - 3 hours, depending on 
the source. This, in turn, suggests that the location of the y-ray 
emission region could be constrained as within the broad-line re- 
gion, thus disfavoring the other possibility of a dissipation zone 
beyond the infrared torus. These upper limits are not conclusive 
yet, because it is still possible to invoke exceptionally great val- 
ues of the Doppler factor to reconcile the observations with the 
hypothesis of the dissipation zone being located at parsec scale. 

Incidentally, we have noted that the shortest upper limit mea- 
sured in this work (t < 2.3 hours) can be set in the case 
of PKS B1222-1-216, which was also recently detected above 
100 GeV by MAGIC (70-400 GeV, Mariotti et al. 2010, Aleksic 
et al. 201 1) and FermifLKT (Neronov et al. 2010), the latter by 
integrating the data collected over several days. The GeV flux 
resulted to be also extremely variable, with doubling time scale 
of the order of 10 minutes, which is now the shortest time scale 
ever detected in a FSRQ (Aleksic et al. 2011). The upper limit 
estimated in the present work is in agreement with the findings 
of the MAGIC Collaboration, although the latter work refers to 
an observation performed on 2010 June 17, at the end of the 
time period analyzed in the present work (see Fig.|2] left panel), 
while the shortest variability measured by Fermi/LAil reported 



here refers to the end of 2010 April, at the beginning of the light 
curve shown in Fig.|2] 

The detection at hundreds of GeV, together with fast flux 
variability, poses new problems. As known, the BLR in FSRQs 
is a rich environment of soft photons, which in turn can severely 
limit the escape of y rays with energies above tens of GeV be- 
cause of pair production (e.g. Liu & Bai 2006). Therefore, these 
photons of hundreds of GeV should come from zones outside 
the BLR, possibly around the molecular torus as suggested by 
Tanaka et al. (2011) on the basis of the analysis of y-ray spec- 
tra, but at such distances from the central black hole, it is ex- 
pected that the blob - if linearly expanding with a constant ratio 
r/R ~ 0.1 - 0.25 - has reached a size so large to exclude the 
possibility of changes over hours time scales or less. A possible 
solution has been recently suggested by Tavecchio et al. (201 1), 
who proposed to explain these observation with something sim- 
ilar to the structured jet proposed for BL Lacs. 
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